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Algorithms are developed for generating a class of exact braneworld cosmologies, where a self- 
interacting scalar field is confined to a positive-tension brane embedded in a bulk containing a 
negative cosmological constant, ft is assumed that the five-dimensional Planck scale exceeds the 
brane tension but is smaller than the four-dimensional Planck mass, ft is shown that the field 
equations can be expressed as a first-order system. A number of solutions to the equations of 
motion are found. The potential resulting in the perfect fluid model is identified. 
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I. INTRODUCTION 
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The possibility that our observable universe may 
be viewed as a domain wall embedded in a higher- 
dimensional space has received considerable attention re- 
cently. Motivated by developments in superstring and 
M-theory B, it is assumed in this scenario that the 
standard model interactions are confined to a (3 + 1) 
dimensional hypersurface, but that gravity may propa- 
gate through the 'bulk' dimensions perpendicular to the 
brane. This change in viewpoint has important conse- 
quences for early universe cosmology and, in particular, 
for the inflationary paradigm. 

Various cosmological aspects of branes embedded in 
five dimensions have been investigated | The effects 
of including a scalar field in the bulk have also been stud- 
ied [p^J15[. In this paper, we consider a co-dimension 
1 brane with positive tension, A, embedded in vacuum 
Einstein gravity with a negative cosmological constant. 
This corresponds to the scenario introduced by Randall 
and Sundrum where the extra dimension is infinite 
We focus on the region of parameter space defined by 
A 1 / 4 <C TO5 <C m.4, where m^^ represent the four- and 
five-dimensional Planck scales, respectively. 

One approach to determining the expansion of the 
brane is to project the five-dimensional metric onto the 
brane world- volume {7|,[L5|. In effect, this is equivalent to 
solving Einstein's field equations, together with appropri- 
ate jump conditions, for a negative cosmological constant 
and an energy-momentum tensor restricted to the brane 
j| . The field equations admit a first integral that may be 
interpreted as a generalized Friedmann equation 
We emphasize, however, that the dynamics could differ in 
a compactified scenario, due to the effects of radion stabi- 
lization. If the cosmological constant and brane tension 
are related in an appropriate fashion, this equation takes 
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where p represents the energy density of matter on the 
brane, e is a constant and we have assumed that the four- 
dimensional cosmological constant is zero. The last term 
on the right-hand side of Eq. ( |1 . 1| ) behaves as 'dark ra- 
diation' and arises due to the backreaction of the bulk 
gravitational degrees of freedom on the brane Pj5|,p|,^0[ . 
At sufficiently low energies, p <C A, the standard cosmic 
behaviour is recovered and the primordial nucleosynthe- 
sis constraint is satisfied provided that A > (lMeV) . 

Quantum gravitational effects become important if the 
energy density on the brane exceeds the five-dimensional 
Planck scale. In this case, the assumption that matter is 
confined to the brane may become unreliable. However, 
if A <C m|, there is a region of parameter space, corre- 
sponding to A <C p <C mf, where the classical solution 
is still valid, but where the quadratic correction becomes 
significant. In particular, this term can play an impor- 
tant role during inflation |p|,pT[. Cosmological inflation 
has played a central role in studies of the very early uni- 
verse. (For recent reviews, see, e.g., Refs. [|Po||l7]]). It is 
therefore important to study the inflationary dynamics 
of the braneworld scenario. In the standard inflationary 
cosmology, the universe is dominated by a scalar 'infla- 
ton' field self-interacting through a potential, V((f>), with 
an energy density p = 4> 2 /2 + V. Inflation proceeds if the 
potential energy of the field dominates its kinetic energy. 
Maartens et al. have recently considered the case where a 



*In this paper we assume that the world-volume metric of 
the three-brane is the spatially flat, Friedmann-Robertson- 
Walker (FRW) line element with scale factor, a(t), and Hub- 
ble parameter, H = a/a. A dot and prime denote differenti- 
ation with respect to cosmic time, t, and the scalar field, <j>, 
respectively. 



1 



single inflaton field is confined to the brane and have de- 
rived the necessary criteria for successful inflation when 
the slow-roll approximation is valid ||. I n ge neral, the 
quadratic term in the Friedmann equation (1.1) results in 
an enhanced friction on the field, implying that inflation 
is possible for a wider region of parameter space than in 
the standard cosmology (||],[l(],|ll| . 

When the scalar field is confined to the brane, energy- 
momentum conservation implies that its equation of mo- 
tion has the standard form: 



ij> + 3H<j) + V'((t)) = 0. 
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Eq. (1.2) can be expressed in terms of the energy density 
such that 
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and it follows that a given inflationary braneworld model 
is specified by a sol ution to the Friedmann-scalar field 
equations (1.1) and (1.3). 

In general, however, it is very difficult to solve this 
system of equations, even when the standard slow-roll 
assumptions, <p 2 <C V and \<f)\ <C H\<j)\, are made. It 
is therefore important to develop ge nera ting te chn iques 
for finding exact solutions to Eqs. (1.1) and (1.2) and 
this is the purpose of the present work. During inflation, 
the dark radiation redshifts rapidly and soon becomes 
dynamically negligible and we therefore consider models 
where e = 0. On the other hand, we do not assume that 
the slow-roll approximation is necessarily valid. Thus, 
such a study has direct applications to the final stages 
of inflation, where the kinetic energy of the inflaton field 
inevitably becomes significant. Depending on the form 
of the potential, this may happen w hen the quadratic 
correction to the Friedmann equation (1.1) is still impor- 
tant. The search for exact scalar field braneworlds is also 
important because it allows one to classify the possible 
types of behaviour that arise in these universes and to 
uncover the generic characteristics of such models. 



II. FIRST-ORDER FIELD EQUATIONS 

From a particle physics perspective, it is natural to 
begin by specifying the functional form of the poten- 
tial. However, even for simple choices, such as expo- 
nential or power law potentials, analytical progress can 
not be made. An alternative route, originally introduced 
within the context of standard chaotic inflation [ pH , is 
to specify the time-dependence of the scale factor, a(t). 
In the braneworld scenario, given a(t), one may deduce 
p(t) from Eq. (1.1) and hence <p(t) from Eq. ( |l.3[ ). Inte- 
grating yields (j>(t). The time-dependence of the poten- 
tial follows immediately from the definition of the energy 
density and the form of V((f>) is deduced by inverting <p{t). 
The drawback of this approach is that such an inversion 
is not always possible and moreover a realistic potential 
does not necessarily result. 



In principle, these problems are partially avoided by 
first specifying an invertible form for 4>(t) |19[, although 
we do not explore this possibility further here. Our ap- 
proach is to note that during inflation, the scalar field 
rolls monotonically down its potential. Thus, Eq. (1.3) 
may be expressed in the form 
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The formal limit, A — > oo, corresponds to the standard in- 
flationary scenario. In this case, Eqs. (1.1) and (2.1) can 
be written in the particularly simple 'Hamilton-Jacobi' 
form 
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where the Hubble parameter is vie wed implicitly as a 
function of the scalar field. Eq. (2.2) may be integrated 
to yield the scale factor: 



a(4>) = exp 
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and this implies that the cosmological dynamics is de- 
termined, up to a single quadrature, once the functional 
form of H (<fi) has been specified. It has been suggested 
that H((j>) should be viewed as the solution generating 
function when analysing inflationary cosmologies p2[ . 
The advantage of such an approach is that the form of the 
potential is readily deduced from Eq. fl2.4p . The ques- 
tion that naturally arises within the braneworld context, 
therefore, is whether there exists an analogous generating 
function when the quadratic correction in the Friedmann 
equation (p~l) is relevant. 

To proceed, we define a new function y(4>): 
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where the restriction y 2 < 1 must be imposed for the 
weak energy condition to be satisfi ed. Substituting Eq. 
( |2.6|) into the Friedmann equation (1.1) implies that the 
Hubble parameter is given by 
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The right-hand side of Eq. (2.7) can be expanded as a 
geometric progression: 
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implying that y is proportional to the Hubble parameter 
in the low-energy limit , y — > (p/A — > 0). 

Substitution of Eq. (2.6) into the scalar field equation 



(2.1) implies that 
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and t he p otenti al is given in terms of y((j>) by combining 
Eqs. and (p7o|) : 



and the potential takes the simple form 
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The time-dependence of the scalar field is determined by 
the integral 
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It follows from Eqs. (2.7) and ( |2.9| ) that the scale factor 
satisfies 
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Eq. ( 2.11 ) is formally identical to Eq. (2.2) and may be 
integrated to yield the scale factor in terms of a single 
quadrature with respect to the scalar field: 



a{4>) = exp 
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Finally, the dependence of the scalar field on cosmic time 
is deduced by evaluating the integral 
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and inverting the result, where to is an arbitrary integra- 
tion constant. 

Thus, we have re duce d the second-order system of 
equatio ns (1.1 ) a nd ( 1.2) to the non-linear, first-order 
system (2.9) and (2.11) by employing the scalar field as 
dynamical variable. 

Further insight may be gained by defining a second 
new function b(<f>): 



y — tanhfe. 
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It follows immediately from Eqs. ( f2.7| ) and (2.9) that the 
Hubble parameter is given by 



H = ( ^) 1/2 sinh2b 
\3m 2 J 

and that the scalar field varies as 
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Eq. ( 2.16 ) is formally equivalent to Eq. (2.3). Hence, 
the scale factor is given by 
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The function b plays the equivalent role to t hat of the 
Hubble parameter in the field equation ( |2.3| ). Indeed, 
b((f>) o c H{4>) in the low-energy limit. It follows from Eq. 
(§T6|) that b(t) is a monotonically decreasing function of 
cosmic time. 

To summarize this Section, we have found that the 
braneworld field equations can be rewritten after suitable 
redefinitions in a way that directly extends the Hamilton- 
Jacobi form of the standard scalar field cosmology. How- 
ever, the physical interpretation of the variables is differ- 
ent in the two cases. Nevertheless, this correspondence 
implies that similar techniques may be employed to find 
exact braneworld cosmologies. In particular, when the 
functional form of the parameter y(cj>) is known, the po- 
tential is determined in terms of this function and its first 



derivative. The function a((f>) foll ows f rom Eq. ( 2.12 ) and 
4>(t) follows by evaluating Eq. (2.13) and inverting the 
result. An alternative method for solving the field equa- 
tions is to specify the for m of b(4>) and t o then evaluate 
the two integrals in Eqs. ( |2.17| ) and ( 2.1S ). This is equiv- 
alent to determining the kinetic energy of the scalar field 
as a function of the field itself. A related technique was 
recently employed in the standard inflationary scenario 
p3| . The ad vanta ge of this approach is that the inte- 
grand in Eq. ( 2.17 ) can be expanded as a power series in 
b. 



III. EXACT BRANEWORLDS 

We now find exact braneworld models by employing 
the above techniques. Firstly, we consider the ansatz 
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where C is an arbitrary constant. Substitution of Eq. 
( |3.lD into Eq. ( (2.10 ) implies that the potential is given 

by 



(3.2) 



V = ±(6-C a ) cosech 2 

3 V ' \ 7714 



Thus, we require C 2 < 6 for the potential to be positive- 
definite and we consider this region of parameter space 
in what follows. Evaluating Eq. ( 2.13| ) implies that 



3 



= {^x) c^ cosh 



t - t. 



2irC 



m 4 



(3.3) 



and substituting Eqs. (3.1) and (3.3) into Eq. ( 2.1 2| ) 
implies that the scale factor is given by 
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Without loss of generality, we may choose to = 
— [3ml/ (AttXC 4 )] 1 ' 2 such that the origin of time corre- 
sponds to a vanishing scale factor. 

It can be verified by direct substitution that the solu- 
tion (3.4) is a scaling solution, in the sense that the ki- 
netic and potential energies of the scalar field redshift at 
the same rate as the brane expands. The field behaves as 
a perfect fluid with an effective equation of state, p = up, 
where the barotropic index is given by 
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Thus, we have found the scalar field model equivalent to 
the perfect fluid cosmology presented in Ref. Q|. The 
solution is interesting because it reduces to the power- 
law cosmology driven by an exponential potential in the 
low-energy limit. Inflation proceeds indefinitely into the 
future if C 2 < 2 and the expansion decelerates for C 2 > 
2. However, at early times, the asymptotic behaviour is 
a cx t x / c , and inflation proceeds for a finite time for 
C 2 < 1. In this limit, the potential is of the form V cx 
(j>~ 2 , where the constant of proportionality determines 
the power of the expansion. 

At late times the function y given in Eq. (|3.l| ) asymp- 
totes to an exponential form. It is therefore of interest 
to consider a second ansatz 
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The potential is negative for y 2 > (6 — C 2 )/6, has a single 
maximum located at y 2 = (6 — C 2 )/(6 + C 2 ) and expo- 
nentially decays to zero from above as y 2 — > 0. Although 
the potential is negative to the right of the maximum, 
the solution exists for all <f> > because the initial mag- 
nitude of the field's kinetic energy is sufficiently large for 
it to move over the maximum and reach <j) — ^ +oo. 

We now consider the formulation of the cosmological 
brane equations summarized in Eqs. ( |2.14 )-(2.19). One 
ansatz that can be invoked is 
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where A is a constant. From Eq. ( [2.16 ), this represents 
a model where the kinetic energy of the scalar field is a 
constant for all time: 



4> = 0o - A (t - t Q ) . 
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The scale factor an d po tentia l of the field are readily 
deduced from Eqs. (2.17) and (2. IS), respectively: 
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Early times correspond to the limit <fi 3> X^^m^/A, 
where the potential has an asymptotically exponential 
form and inflation may proceed for a wide range of pa- 
rameter space. Inflation ends within a finite time, how- 
ever, and the field eventually falls into its minimum at 
V m i n = — A 2 / 2. The expansion of the universe is then 
reversed and the subsequent collapse enables the field to 
move up the other side of its potential. 
Another solvable model is defined by 



that is valid for all time, where we assume implicitly that 
<j) > and that C 2 < 6. Integrating Eqs. (^ and fl2.12| ) 
implies that 
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and the origin of time corresponds to T = 1. The po- 
tential of the s calar field is deduced by substituting Eq. 
@ into Eq. ( |2~10| ): 
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b = p 2 <p 2 
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where pi is an arbi trary constant. Substituting Eq. 
( |3.15| ) into Eq. ( 2.1g ) implies that the potential is given 

by 
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This potential has a single maximum at V(<fi = 0) = 
and two minima at 
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Integration of Eq. ( [2.16] ) implies that the time- 
dependence of the scalar field is given by 



t> exp 



\Xp\m\ x 1/2 
3tt 



(t - t ) 



(3.18) 



4 



Finally, the growth in the scale facto r can be expres sed a s 
a power series by substituting Eq. ( 3.15 ) into Eq. ( 2.17 ) 
and integrating: 
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Before concluding, we develop in the next Section an 
algorithm that generates a new braneworld cosmology 
from a known solution such as those presented above. 



IV. A THIRD ALGORITHM 

Sin ce bit) is a monotonically decreasing function, Eq. 
( 2.16 ) may be rewritten as 



2\ 1/2 



b. 



Defining a new function 



c = exp 



and a new time parameter 
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implies that Eq. (4.1) may be then expressed in the form 
of a one-dimensional Helmholtz equation: 
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where the effective potential is uniquely determined by 
the kinetic energy of the scalar field: 
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The importance of Eq. (4.4) is that if a particular solu- 



tion, 01(77), is known, the general solution can be written 
in terms of a single quadrature with respect to this solu- 
tion M: 
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where k is an arbitrary constant. It is this feature that 
form the basis of the algorithm. Suppose a particular so- 
lution to the braneworld field equations has already been 
found, i.e., that {b(t) , <fi(t)} are known. The func tion, 
01(77), can in principle be evaluated from Eqs. (4.2) and 
(4.3). Eq. ( |4.6| ) then yields the new solution for 0(77) and, 
hence, 6(77) from the definition (|4.2|). The form of 4>{rj) is 



identical in both solutions, but the potential of the scalar 
field is different in the new solution. It is given by 
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and 



The form of V(<p) follows from Eqs. ( |4~7| ) and ([4~ 
the scale factor is deduced as before from Eq. ( 2.17 ). 

In effect, each particular braneworld cosmology is 
twinned with a second solution. In practice, it ma y not 
always be possible to evaluate the integrals (|4.2j) and 



(4.3), but it would be interesting to explore models of 
this nature further. We also remark that in the low- 
energy limit, c varies as a power of the scale factor of 
the brane. Thus, the above discussion is also relevant to 
standard scalar field cosmology. 



V. CONCLUSIONS 

In conclusion, we have presented algorithms for solv- 
ing the braneworld Friedmann equation for a single, self- 
interacting scalar field confined to a brane with positive 
tensione3 embedded in five-dimensional, vacuum Ein- 
stein gravity. The formalism does not assume the slow- 
roll approximation and is valid during inflation when the 
field is monotonically rolling down its potential. It does 
not apply if the field is oscillating about a minimum. 
A number of new exact solutions were found, includ- 
ing the scalar field model that corresponds to a perfect 
fluid. The algorithms may be viewed as generalizations 
of the Hamilton-Jacobi formalism that has played a cen- 
tral role in analyses of the standard chaotic inflationary 
scenario |2(J. In the latter case, the Hamilton-Jacobi 
formalism provides the necessary framework for estab- 
lishing the precise correspondence between the potential 
of the scalar field and the scalar and tensor perturba- 
tion spectra that are generated during inflation |l7j]. It 
would be interesting to employ the techniques developed 
above to establish the equivalent correspondence in the 
braneworld scenario. 
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